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_ ■ Abstract. We study the evolution of the first stars in the universe (Population III) from the 

early pre-Main Sequence (MS) until the end of helium burning in the presence of WIMP dark 

. matter annihilation inside the stellar structure. The two different mechanisms that can provide 

' this energy source are the contemporary contraction of baryons and dark matter, and the capture 

\l \ of WIMPs by scattering off the gas with subsequent accumulation inside the star. We find that 

' . the first mechanism can generate an equilibrium phase, previously known as a dark star, which is 

' transient and present in the very early stages of pre-MS evolution. The mechanism of scattering 

and capture acts later, and can support the star virtually forever, depending on environmental 

<^ |i characteristic of the dark matter halo and on the specific WIMP model. 
I ■ 
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1. Introduction 

Within the scenario of a ACDM cosmology, it has been recently recognized that if 
the dark matter (DM) dominant component is a weakly interacting massive particle 
(WIMP), its annihilation could play a relevant role on the formation and evolution of 
the first baryonic objects in our universe. Spolyar, Freese, and Gondolo (2008) noticed 
fSl that during the proto-stellar phase, the cooling of the baryonic gas could be overcome 
, by the energy deposition following the annihilation of DM concentrated in the star for- 
mation site. This is a consequence of the peculiar formation characteristics of the first 
OO stars, at the center of a minihalo whose gas cooling is dominated by the little efficient 
primordial chemistry; the authors suggested that this phase (called a dark star) could 
^ , prevent the formation of the first stars, and be a new phase of stellar evolution. locco 
(2008) and Freese, Spolyar, and Aguirre (2008) noticed that, if a star does eventually 
form, the process of WIMP capture by scattering could be so efficient that their subse- 
ts quent accumulation and annihilation inside the celestial object may provide an energy 
source comparable or even exceeding its nuclear luminosity. Motivated by these works, 
locco et al. (2008) and Freese et al. (2008b, c) have studied the early pre-MS phase, in 
which the baryonic structure is sustained by the annihilation of the DM accreted inside 
it by gravitational contraction, with the help of numerical codes. Both the groups find 
this phase is transient, although the techniques adopted are different and the details of 
the treatment lead to different duration estimates; they both conclude, however, that the 
collapse must continue at the end of this process, which for the sake of simplicity we call 
the Adiabatic Contraction {AC) phase. In locco et al. (2008) (hereafter, 108), we have 
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also studied the pre-MS phase of these stars in presence of annihilating scattered and 
captured DM {SC phase) and followed the evolution of stellar models of different mass 
until the end of the helium burning, in different DM environments. Wc find the duration 
of the MS is dramatically prolonged, up to a potentially everlasting phase (depending 
on the choice of parameters, see later): the energy released by the annihilating DM can 
support the core at temperatures low enough that nuclear reactions arc never ignited. 
Yoon, locco, and Akiyama (2008) and Taoso et al. (2008) also studied the SC phase, 
confirming our results with different codes and carrying their analysis further. These 
proceedings arc based on the results obtained by 108, to which we address the reader 
for detailed referencing and more quantitative details: here we aim to a more qualitative 
description of the physical processes at the basis of this class of objects. 



2. Adiabatic Contraction phase 

Early stars are thought to form in halos of M/j r^lO^ Mq and virial temperature T<10^ 
K at rcdshift Z'^20. The primordial, metal free composition of the gas, the absence of 
strong magnetic fields and turbulence make early star formation very different from the 
one in the older universe. Simulations tell us that the result of such peculiar environment 
are massive (30-300 Mq) stars that form in the very center of the halo (accurate review 
and referencing in the First Stars 111 proceedings, 2008). The collapse of the baryonic 
material "pulls" also the coUisional DM towards the center, thiis contributing to the 
build-up of a central "spike" of DM; if that is indeed made of self-annihilating particles 
(as WIMP neutralinos are), its higher concentration in the center of the halo causes a 
huge enhancement of the annihilation term: 

dV m,j^ ' 

with p the local DM density, the neutralino mass which wc take to be lOOGeV-, 
{av) the thermally-averaged annihilation rate. It is worth noting that this energy term 
depends on the self-annihilation rate, whose value can be quite safely established through 
cosmological arguments as {av) = 3 x 10~^^ cm'^ s^^, see for instance a recent review on 
particle DM (e.g. Bertone et al. 2005). A fraction (1 — /) of the whole energy is emitted in 
the form of neutrinos, and therefore lost by the system, whereas all the other products of 
annihilation can quickly thermalize inside the gas already at densities n ^IQ^'^ jj= j cva'^ , as 
shown by Spolyar et al (2008); we take / ~2/3, a typical value for a neutralino annihila- 
tion. In order to study this peculiar object we have modified the Padova Stellar Evolution 
code to take into account the energy released by DM annihilation. The initial models 
have been prepared as follows: (i) the "stellar" profile (namely the baryonic structure) 
has been obtained by "pumping" a Zero Age Main Sequence star with an artificial energy 
source (as much as permitted by the stability of our numerical code); for a lOOM© star, 
the resulting structure is an object of radius R^,=1.2x lO^^'cm and effective temperature 
T ~5xl0^ K; (m) the DM profile has been obtained in the approximation of an adia- 
batically contracted DM profile (from an original NEW) which has been matched to the 
baryonic structure in the center, whic;h dic;tates the gravitational potential. For details 
on the characteristics of the halo and on the adiabatic contraction approximation used 
we address the reader to 108. This adiabatically contracted DM profile represents the 
initial model, at the time we start our analysis: we have implemented in our code a rou- 
tine which allows to follow the adiabatic contraction of the DM inside our stellar object, 
account for the energy released by its annihilation, and so self-consistently describe the 
evolution of the DM coupled to the baryons. In this way we can follow the contraction 
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Figure 1. In the left panel: DM density profiles, truncated at the stellar radius, for a IOOM0 at 
different times during the AC phase. In the right panel: evolution of the effective temperature 
for all stellar models, benchmarked by the relative contribution of DM annihilation to the total 
luminosity of the star. See text for details 

of the protostar toward the Main Sequence including all relevant energy sources (gravi- 
tational, DM annihilation, and possibly nuclear) and adopt a complete equation of state 
for the gas. 

For the adopted AC parameters, DM annihilation is able to support the baryonic 
structure in the very early phases of gravitational contraction, in all our selected stellar 
models (5-600 M©) The protostars contract "pulling in" more and more DM until an 
equilibrium is reached in very short times, C'(lO^yr). 

However, this phase is doomed to be transient. As the annihilation proceeds, the star 
contracts because of the shrink of DM luminosity caused by its consumption. But this 
further contraction is not able to restore the previous DM density profile, and the star 
finds its equilibrium at a lower luminosity. In this phase the star descends slowly along 
the Hayashi line. Besides consumption due to annihilation, the star continuosly looses the 
external shells of DM, because of the different response af baryons and DM to the grow- 
ing gravitational potential. Eventually the contraction of DM becomes not fast enough 
to restore an "efficient" cusp and the DM luminosity is no more able to balance the 
stellar energy losses. At this point the star keeps on contracting on a typical the Kelvin- 
Helmholtz time-scale, terminating the AC phase and moving toward the main sequence. 
By defining the duration of this phase tac as the the time needed from the DM an- 
nihilation to scale from 100% to 50% of the total luminosity of the object, we find 
typical values of order tag '^lO'^yr, ranging from tac=2.1 x lO'^yr for a GOOM© star up to 
tac=1-8x lOVr for a 9Mq star. In Figure[T]we show the effective temperature evolution 
for different stellar models and the DM profile inside the baryonic structure for a lOOM© 
star at different times during the stellar collapse. The benchmark values correspond to 
the fractional contribution of DM annihilation to the total luminosity. 

3. Scattering and Capture phase 

The other physical mechanism able to concentrate DM inside a star is the capture of 
WIMPs by means of elastic scattering with the gas particles that constitute the object. 
The captured WIMPs thermalize with the gas and eventually reach a thermally relaxed 
state inside the star: the density profile is 

nAR)=n^^eM-RVrl), -g^, (3.1) 
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Figure 2. In the left panel: HR diagram for metal-free stars in presence of a DM environmental 
density p=10'^^GeV, compared to their normal position in the HR; sold curve labeled NO SC. 
For all masses, the ZAMS is never reached, as stars are entirely supported by DM annihilation. 
In the left panel: stars in the white region, above the solid line never ignite nuclear reactions, 
being "frozen" in their evolution. In the shaded area they do evolve at different rates, being only 
partially supported by DM burning. See text for details. 



where and pc are the stellar core temperature and density, respectively, and C* the 
WIMP capture rate, discussed in the following. It can be seen that is in general much 
smaller than the star core, ~10^ cm for a lOOM© star. For each volume element, 
the DM annihilation energy released is given by Equation 12.11 however, capture and 
annihilation reach the equilibrium on timescales much shorter than the stellar lifetime 
(see locco (2008) and 108 for a more detailed discussion of this issue), and the amount of 
energy released inside the star, the "dark luminosity" Ldm reads LDM=fm^Ci,, being 
dictated by the capture rate, C*. The latter has the characteristic of a scattering term, 
and can therefore be cast as: 

C.(xM,t;L — • (3.2) 

with (To the elastic scattering cross section between WIMP and baryons, p the DM density 
outside the star, Vesc the escape velocity at the surface of the star, and the subscript 
* referring to stellar quantities throughout this paper. We refer to 108 and references 
therein for a more quantitative discussion of the capture rate, but we wish here to stress 
on few peculiarities of the SC process. The energy is indeed provided by means of DM 
annihilation in the bosom of the star; however, the bottleneck of the capture/scattering 
mechanism is given by capture: the star cannot burn more WIMPs than it accretes. 
Therefore, at the equilibrium, the annihilation is dictated by the capture and shows the 
parameter dependence of a scattering process. Also, this process is sensitive to the WIMPs 
that stream through the star, and thus the halo plays the role of a "reservoir", virtually 
unexhaustablcQ. It depends on ctq rather than (av), and linearly on the environmental 
DM density p. However, when the object is in the phase, the SC mechanism is playing 
little role because of the small capture rate (the star is an inefficient "net" for capturing 
DM) and of the very long equilibrium timescale at these stages, much longer than the 
Kelvin-Helmholtz time. By contracting, the star enhances the efficiency of the SC process 
and DM annihilation contributes to the stellar luminosity. 

For f7o=10~^^ (at the level of the current upper limit for the spin-dependent elastic 

f The DM mass-energy content of the halo is much bigger than the energy used by the 
star: only approximately lOM0 = lO~^Mh are entirely converted into energy in lOGyr, at a rate 
10^'^erg/s. Only O.IM© of DM can support a IOOMq star for ten times longer than its natural 
lifetime, lO^yr. 
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scattering cross section (Desai et al., 2004) ), and an environment density /5=10^^GeV/cm'^ 
(a likely value achieved around the star, see Figure [1]), the energy provided by the DM 
annihilation inside the star is able to entirely support all our models before they get to 
the ZAMS, therefore not igniting nuclear reactions. Their locus on the HR diagram can 
be observed in Figure [2l compared with their ZAMS one. Lower values of p or ctq make 
the DM contribution smaller, and the star can contract increasingly, thus progressively 
enabling the ignition of nuclear reactions. Stars in the white region of the right panel 
of Figure [1] do never ignite nuclear reactions; stars below the solid curve, in the shaded 
area do. The closer to the curve, the longer their lifetimes will be, as energy contribution 
required by nuclear burning is lower, and their chemical evolution is slower, see 108. 

4. Conclusions 

We have studied the effects of WIMP dark matter annihilation from the late stages of 
stellar formation through the end of helium burning. Addressing the reader to locco et al. 
(2008) for quantitative details, here we wish to stress the fundamental physical difference 
between the SC mechanism, which needs a weak process to mediate the capture and 
thermalization of DM particles (and therefore the DM density profile inside the star), 
and the AC one, where the DM density profile is dictated by gravitational contraction 
only. The difference in the two physical processes dictates diverse characteristics for the 
two phases: the AC is transient, and takes place early in the pre-Main Sequence evolution 
of the star, when it is almost at the proto-stellar stage. The SC phase becomes active 
when the star is at the bottom of the Hayashi track or later in the pre-MS evolution, 
and it can dramatically extend the stellar life, up to orders of magnitude more than its 
standard lifespan, depending on parameters. 

As a note added in proof, we address the reader to a paper very similar to this one, 
Freese et al. (2008d) in these proceedings; and also to the analysis of Natarajan, Tan, 
and O' Shea (2008) who find results in good agreement with Spolyar et al. 's, studying 
gas and dark matter profiles from cosmological simulations of first star formation. 
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